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ABSTRACT 

We present measurements of the cosmic microwave background (CMB) radiation temperature anisotropy in 
the multipole range 2000 < £ < 3000 from the QUaD telescope's second and third observing seasons. After 
masking the brightest point sources our results are consistent with the primary ACDM expectation alone. We 
estimate the contribution of residual (un-masked) radio point sources using a model calibrated to our own 
bright source observations, and a full simulation of the source finding and masking procedure. Including this 
contribution slightly improves the \ 2 - We also fit a standard SZ template to the bandpowers and see no strong 
evidence of an SZ contribution, which is as expected for w 0.8. 
Subject headings: cosmology: cosmic microwave background, cosmology: observations 



1. INTRODUCTION 

Observations of the cosmic microwave background (CMB) 
anisotropy at angular scales of several arcminutes or larger 
(£ < 2000) have been used to constrain param eters of the 
ACDM cosmological mo del to high precision (ICastro et al.1 
2009: lDunklev et al.l2009l) . At these larg er angular scales, the 
anisotropic power is dominated by the primary CMB from the 
surface of last scattering. At smaller angular scales {£ > 2000) 
the primary anisotropy is exponentially suppressed by dif- 
fusion in the primordial plasma and the structure becomes 
dominated by foreground emission and secondary anisotropy 
generated by intervening large scale structure. Measuring 
the secondary anisotropy introduced by the thermal Sunyaev- 
Zel'dovich effect (SZE) has been of particular interest. The 
magnitude of the SZE power is a sensitive and indepen- 
dent probe of the amplitude of density perturbations, scal- 
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ing as erg (Komatsu&Seljak 2002). 

Previous measurements of the small angular scale CMB 
anisotr opy at 30 GHz by CBI dReadhead 1 et alj 120041) and 
BIMA (Dawson et al. 2006) claimed a significant excess over 
the ACDM expectat ion at multipoles of £ > 2000. The 
ACBAR experiment dReichardt et a l. 2009) subsequently re- 
ported a ~ 1-cr excess at 150 GHz at similar scales. At- 
tributing this excess power to the SZE alone implies erg ~ 
1. This value is in conflict with the WMAP 5 -year results 
dDunklev et al J 120091) and recent X-ray me asurements of the 
cluster mass function dVikhlinin et al.l2009l) . which both yield 
values of erg ss 0.8. For the latter value of erg, the SZE power 
at 30 GHz is expected to be comparable to the primary CMB 
at multipoles of £ « 2500 but at 100 and 150 GHz will be sub- 
dominant at multipoles of £ < 3000. The results presented in 
this work cover a multipole range of 2000 < £ < 3000 and are 
the highest sensitivity to date at these scales. 

The QUaD telescope is a millimeter-wavelength bolomet- 
ric polarimeter located at the South Pole. QUaD operated 
during the austral winters of 2005 to 2007. Details of the 
Q UaD instrument , calibrations and performance can b e found 
in lHinderksetail d2009l) and |Q' Sullivan et all d2008l) . In this 
paper we present high-£ TT spectra only. Details of the ob- 
servations, data quality, low level processing, map-making 
and power spectrum estimation plus t he full polar i zation anal- 
ysis for £ < 2000 can be found in iPrvke et all d2009l) and 
iBrown et alJ ((2009) — this paper follows the analysis meth- 
ods described there except where noted. 

2. ANALYSIS 

In order to reduce the bandpower uncertainty at high-^ we 
have adopted an optimal signal to noise Fourier plane weight- 
ing step in the power spectrum estimation. As can be seen in 
Figure 7 of Prv keet al.l J2009). the distribution of noise power 
in the two-dimensional Fourier plane is highly non-uniform 
— the atmospheric noise forms a concentrated band around 
the y-axis. Fourier plane weights are calculated as 



S 2 

J CMB 



(SCMB+AO 2 



(1) 



where Scmb an d N stre the ensemble averages of the signal 
and noise simulation two dimensional auto power spectra. 
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FIG. 1.— TTbandpower values at 100 GHz (top) and 150 GHz (bottom) 
versus ACDM. We show results using the same data from two independent 
pipelines (see text); the points are offset for clarity. The error bars include 
ACDM sample variance and noise only. The triangles indicate the coherent 
shift in the errorbar end points which result when the beam and absolute 
calibration are simultaneously pushed up/down by l-er. 
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FIG. 2. — Jackknife 7Tbandpower values for both the 100 GHz (top) and 
150 GHz (bottom) deck-jackknife maps. These spectra are consistent with 
null and indicate that the signal seen in FigurefT|originates on the sky. 

Since the CMB signal is expected to be uniformly distributed 
in azimuth angle, downweighting localized regions of high 
noise in the Fourier plane will not bias the results so long as 
the weighting is independent of the actual data values. This 
weighting has a dramatic effect on the bandpower uncertainty 
at high-£ — for 150 GHz the error is suppressed by as much 
as an order of magnitude in the range of 2500 < £ < 3000. 

In this paper w e also use two enhancements over the 
iPrvke et all ( 120091) analysis introduced in lBrown etail (120091) . 
Firstly we replace the field differencing operation used to re- 
move ground contamination with a template removal tech- 



nique that doubles our effective sky area, further reducing the 
bandpower uncertainties by a factor of ~ y/2. Secondly the 
beam model used in this paper has been updated to include 
sidelobes measured at the < -20dB level and predicted by 
the ph ysical optics simulations described in |Q' Su llivan et al. 
(2008). The absolute calibration uncertainty is now 7% in 
power. The cosmological model assumed in our simulations 
and analysis is t he WMAP 5-year m odel given in column two 
of Table 2 in iDunkley et al.l (120091) with zero SZE signal, 
hereafter referred to as ACD M. The inpu t simulation spectra 
are generated using CAMB dLewis et al.ll2000h with lensing 
turned on. 



3. RESULTS 

In Figure [T] we present our basic result, the TTband power 
values at 100 and 150 G Hz extending to t = 3000. Pipeline 



A is the pipeline used in IPrvke et al.rd2009[). while Pipe line 
B is an alternate cu rved sky analysis (see lAde et al.ll2008l and 
iBrown et alj|2009b . Both pipe lines are based br oadly on the 
MASTER analysis technique dHivon et al.ll2002l) . The band- 
power uncertainties are calculated from the spread in signal 
plus noise simulation bandpowers assuming the input ACDM 
theory spectrum shown. 

As discussed in detail in Prv ke et aT] d2009h . jackknife maps 
made from differencing independent data sets covering the 
same sky are a powerful test for systematic contamination. In 
Figure|2] we present bandpower values for the deck-jackknife 
— likely our most stringent test. These jackknife spectra are 
consistent with null. 

In the multipole range considered for this analysis, the sky 
power has been suppressed by almost an order of magnitude 
through beam convolution. Thus a small mis-estimate of the 
beam would result in a large, multipole-dependent systematic 
shift in the bandpower values. An under(over)-estimate of the 
beam suppression would result in an under(over)-estimate of 
our bandpowers. The effect of the systematic uncertainties 
on our beam model and calibration is illustrated in Figure Q] 
where we show the result of pushing both up/down simulta- 
neously by l-cr. While systematic uncertainty is significant at 
both frequencies, it is not sufficient to qualitatively change the 
results at 150 GHz. 

Though it is customary to present the bandpower results as 
in Figure[U it is arguably more natural to consider them as we 
do in Figure[3] We calculate \ 2 between the data and the sim- 
ulation distributions using the spread in the signal plus noise 
realizations to construct a bandpow er covariance matrix (see 
IPrvke et al.l2009HBrown et alj2009l) . The resulting x 2 values 
calculated over the I bins presented are shown in Figure [3] 

The bandpowers presented in Figures Q] & [3] were calcu- 
lated after masking bright point sources in the maps as de- 
scribed in Section [4] We detect 7 point sources at > 5(7 
(~ 50mJ) in both the 100 and 150 GHz maps; all of these 
were mat ched with a low-pr obability of chance-association 
to PMN dGregorv et aljl994h or SUMSS dMauch et alj2003l) 
radio sources using NED 1 . The effect of masking them can 
be seen in Figure [3] as the difference between the light and 
dark points. Since radio source populations typically follow a 
power law distribution these sources are only the sparse high 
flux end of an exponentially more numerous low flux popula- 
tion. It is therefore necessary to estimate the residual power 
contribution from the unmasked radio source population. 

1 http://nedwww.ipac.caltech.edu/ 
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FIG. 3.— TT bandpower values at 100 GHz (top) and 150 GHz (bot- 
tom). This Figure shows the same information as Figure ff] but instead of 
using error bars we plot the spread in the signal plus noise simulations versus 
bare points from the data. This is in a sense the fundamental result of the 
analysis — is the data consistent with being a realization of the simulation 
model? The lines show the 16, 50 and 84% points of the simulation distri- 
bution (corresponding to -1, 0, +1 a for a Gaussian distribution). Simulation 
distributions are shown for ACDM alone and for ACDM plus a residual radio 
source foreground (pntsrc). The x values are calculated for the data versus 
the simulation model (see text). We also plot bandpowers calculated without 
masking sources in the maps as light points. 

4. POINT SOURCE SIMULATIONS 

The effect of residual point source contamination in our 
spectra will manifest as both an increase in the total power at a 
given I and an increase in the bandpower fluctuation. Though 
it would be straightforward to estimate the mean power con- 
tribution from a given point source model and flux cut analyt- 
ically the subtle effects of source identification and masking 
in our pipeline would be difficult to account for accurately. 
Moreover, the fluctuations are potentially non-Gaussian. In- 
stead, we explicitly simulate the source population in our 
maps. 

At present, the statistical properties of the sub-Jansky radio 
source population are not well known at 100 and 150 GHz. 
The most useful published data is the W-ban d (94GHz) 
WMAP point source catalog (IWright et al.ll2009l) . However, 
this catalog is only complete at fluxes at or above several Jan- 
sky whereas even the brightest sources in the QUaD maps 
are sub-Jansky. Therefore to predict the radio source distri- 
bution below our detection threshold we use the lde Zotti et al.1 
(2005) extragalactic radio source model. Although this is a 
carefully constructed model using detailed astrophysics, it's 
calibration at lower frequencies makes it a distant extrapola- 
tion at 150 GHz. We rescale the model by 0.7 and 0.6 at 100 
and 150 GHz respectively to match the number of sources ob- 
served in the QUaD maps (see Figure|4]i. Doing so also brings 
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FIG. 4. — The simulation input point source model (solid lines) is derived 
by applying a s imple rescale factor to the de Zotti radio point source model 
(dashed lines) I de Zotti et al. 2005) at each frequency. The rescale factor 
is determined from fits to dN/dS points derived from sources detected in 
our maps (closed circles). At 100 GHz, the re-scaling also brings the de 
Zotti model into better agreement with WMAP W-band dN/dS points (open 
circles). 

the model into better agreement with crude dN/dS points de- 
termined from the WMAP catalog. We have ignored the pos- 
sible contribution from IR (dusty) sources. 

For point sources, unlike the CMB, we do not simulate the 
sky signal at the timestream level. Instead we directly in- 
ject source populations into existing signal plus noise maps 
as beam sized blips. The fluxes are drawn from the source 
population model and the positions are uniformly distributed 
across the map area, i.e. no clustering is assumed. 

While sources of low to moderate flux are abundant those 
of high flux are rare. Leaving the brightest sources in the 
maps leads to huge contamination as we see comparing the 
black and grey points in Figure[3] Thus we must mask sources 
which are detected at high significance. Given a map contain- 
ing three components — CMB, noise and point sources — 
we need an automated and unbiased method for identifying 
the point sources to mask. 

To separate out the point sources in our maps we adopt a 
two-dimensional Fourier space optimal filter (Weiner filter) 
given by 

W = -_ ^ (2) 



Spnt + ScMB + N 

where S pnt is the point source signal. 



We fix S pM to be a 
beam suppressed white-noise power spectrum — represent- 
ing a uniform distribution of point sources — with an am- 
plitude at I ~ 2500 roughly equal to our bandpower value at 
that multipole. The resulting filters are broadly azimuthally 
symmetric in form. They are non-zero in the multipole range 
where we are most sensitive to point source power — i.e. zero 
at the lowest multipoles where CMB signal is dominant then 
rising to a peak near I ~ 2000 for 100 GHz and i ~ 2500 for 
150 GHz before decaying back to zero at the highest multi- 
poles where instrumental noise is dominant. 
The filtered map is then given by 



/ / hi 
= FT I W ■ FT I 



v pvxi 



(3) 



where FT is the Fourier transform operation, m is the original 
map, W is the Weiner filter defined in equation [2] and v pix is 
the pixel variance map — used here to appodize m — esti- 
mated from the timestream rms at the map making stage (see 
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iPrvke et alj|2009b . For source identification we construct sig- 
nal to noise maps as 

' -V 2 fA\ 

s = mv pi ' x . (4) 

Though the filtering operation changes the noise amplitude, 
Vpix provides information about the spatial distribution of the 

noise. We rescale the amplitude of s so that the 16 t ^ 1 per- 
centile point of the pixel distribution is equal to -1. To iden- 
tify sources in s, we subtract a source template — constructed 
from the back-transform of W — from the brightest pixels in s 
and iterate down to a threshold value of 5. This procedure re- 
sults in a source free s-map with pixel values that are close to 
Gaussian distributed in the range ±5cr and a catalog of > 5<r 
sources, where a w 10 mJy at both frequencies. 

We generate a source catalog for both the real maps and 
for each of our point source injected signal plus noise sim- 
ulated maps and calculate the power spectra masking out 
the > 5(7 sources. The resulting distribution of the radio 
source injected simulations, as compared to ACDM alone, 
is plotted in Figure [3] Comparing the data against this new 
model we find that the addition of residual radio sources 
marginally improves the x 2 at both 100 and 150 GHz. The 
bandpowers values, together with their covariance matrices 
and window functions, are available in numerical form at 
http://quad.uchicago.edu/quad 

5. IMPLICATIONS FOR SZE FOREGROUNDS 

Figure [5] shows the QUaD High-£ TT results along with 
measurements from other recent experiments. There is broad 
agreement between the various datasets — amongst each 
other and with ACDM — exc ept at the highest multipoles. 
Here CBI ( Sievers et al. 2009]) clai ms a significant exce ss, 
which was "confir med" by ACBAR dReichardt et al.ll2009t) at 
1-cr, whereas SZA (Shar p et al.ll2009i) is consistent with zero 
power. This intimation of excess power has been attributed to 
SZE signal and used to derive corresponding constraints on 
the amplitude of density perturbations erg. Depending on the 
template used, the CBI (and to a lesser extent ACBAR) data 
imply values in the range 0.9 < cr% < 1 .0, a departure from the 
conventional value of 0.8. 

Foll owing suit, we add the standard Komatsu & Seliak 
(2002) (KS) template — scaled by a single parameter Asz — 
to ACDM plus the residual radio source contribution (Sec- 
tion|4]i and fit this model to the QUaD bandpowers. The scale 
parameter can be related to a value of a% as 

oF«0.8A^. (5) 

The 100 and 150 GHz data are considered both indepen- 
dently and simultaneously. We use the bandpower covariance 
matrix from simulations and add a beam plus absolute cali- 
bration systematic term calculated as 

M' bb , = M w + a 2 (S h C b ) (S V C V ) , (6) 

where Q, is the real bandpower value, a is the absolute calibra- 
tion uncertainty and Sb is the beam uncertainty. We calculate 

ax 2 for A sz as 

x i=(C b -(Cn-AszCD M£J (6 v -(C>* c )-A sz Cf?), 

(7) 

where Q,, Cf c and Cf s are the observed, simulation and KS- 
expectation bandpowers respectively. The sample variance 



component of the covariance matrix is scaled in proportion to 
the level of model signal akin to role of the more conventional 
offset-lognormal transformation. The resulting likelihood is 
still very nearly Gaussian. 

We note that the procedure adopted here neglects the non- 
Gaussianity of the SZE signal and its fluctuation and to cor- 
rectly include this would require that we inject simulated SZE 
sky into our maps — we have not done this. 

We quote the maximum-likelihood value, with the 1-cr un- 
certainties corresponding to the likelihood falloff that encom- 
passes 68% of the total; 95% upper limits are evaluated by 
integrating the likelihood over postive values of Asz- The re- 
sults for 100 and 150 GHz are A sz = 0.9 ± 1 .4 and 1 .0 ± 1 .5 
or Ajz < 3.6 and 3.8 respectively with simultaneous fit values 
of A sz = 1.2±1.2orA sz < 3.3. 

These results are consistent with SZE power at t he expected 
level f or ag = 0.8 but inconsistent with those of ISievers et al.l 
(2009), preferring lower value s of Asz- There is g ood agree- 
ment with the conclusions of ISharp et al.1 (120091) who also 
make a strong argument that the CBI point source estimate 
is in fact too small. While the frequency dependence of the 
SZE makes measurements at 100 and 150 GHz intrinsically 
less sensitive than those at 30 GHz, they are also less prone to 
contamination by radio sources. 



6. CONCLUSIONS 

We have extended the range of the 7T-bandpowers from 
QUaD to I = 3000 using signal to noise weighting to down- 
weight noisy regions of the two-dimensional Fourier plane, 
an improved method for the removal of ground contamina- 
tion and higher accuracy beam modeling. After masking point 
sources detected at high significance in our maps, the spectra 
are consistent with the ACDM expectation alone. We have 
estimated the residual radio source contribution using a phys- 
ically motivated radio source model scaled to fit our bright 
source counts and find the contribution to be small at 100 GHz 
and negligible at 150GHz. 

A small SZE contribution is expected at £ < 3000 for erg = 
0.8 (see Figure |5}. Fitting a standard SZE template spectrum 
to our data results in a best-fit amplitude consistent with the 
expectation. 
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FIG 5.— The QUaD high-l TT results for 2000 < I < 3000 compared against recent results from ACBAR ( Reichardt et af, 2009 |) CBHSievers et alJ2009l) and 
SZA I Sharp et al. 2009) — sp anning the spectral range 30, 100 and 1500 GHz (red, green and blue respectively) — plus WMAP iHinshaw et al. 2009) (black) 
and QUaD iBrown et al. 2009) for i < 2000 (grey). For QUaD, SZA and CBI the estimated residual radio source contribution h as been subtr acted. Some points 
have been slightly offset in multipole for clarity. The data are plotted against ACDM alone and ACDM plus the standard Komatsu & Seliak 12002) template 
assuming two values of erg scaled to all three frequencies. The QUaD data favors the erg =0.8 model with a best-fit scaling of A$z = 1.2 ± 1.2 (see text). 
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